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Abstract. I study the importance of the accurate calibration of photometric 
systems in order to produce meaningful comparisons between the observed colors 
+ magnitudes and model SEDs. Possible sources of errors are discussed and two 
examples are analyzed. I show that well-calibrated Tycho-2 photometry is stable 
and precise enough for such comparisons. On the contrary, the available calibra- 
tions for Johnson UBV photometry yield relative large systematic errors, which 
has prompted me to develop a new, more precise calibration. The advantages 
of multicolor photometry over the standard single-color + magnitude diagrams 
for the derivation of physical properties of stars (elimination of degeneracies, 
inclusion of multiple parameters, avoidance of linearizing approximations, pos- 
sibility of a more precise treatment of errors) are discussed through the use of 
CHORIZOS, a code developed specifically for this purpose. 



1. A brief introduction to synthetic photometry 

The most common way of studying the properties of stellar populations is by 
comparing measured single-color + magnitude diagrams (SCMDs) with the syn- 
thetic magnitudes derived from spectral energy distributions (SEDs) by means 
of a synthetic photometry code. The formula used to compute the magnitude 
for a photon-counting detector is: 

and the equivalent formula for an energy-integrating detector is: 

The quantities in those formulae are: 

• fx (A) is the SED of the object. 

• /A,ref(A) is the SED of the reference spectrum. The most common choices 
are Vega, a constant in fx , and a constant in f v (see Fig. ^) . 

• -P(A) is the total-system sensitivity curve. 
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• ZPp (or ZPp) is the zero point for filter P. 
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Figure 1. Reference sp ectra for magnitude systems: the Vega spectrum of 
Bohlin & Gill ilandl l)2004|) . the constant f\ spectrum for the ST system, and 
the constant /„ spectrum for the AB system. The latter two are normalized 
to have the same flux (in photons s _1 cm" 2 A" 1 ) at A = 5493 A (dotted line). 
The non-normalized photon-counting sensitivity curves for the Tycho-2 Bt 
and Vt filters are also plotted. 

A magnitude system for a filter P is defined by giving both f\ tTe {(X) and 
ZPp. Common choices are: 

• The VEGAMAG system uses Veg reference spectrum and has ZPp 
= by definition JSTScJIlflQfl . 

• The STMAG system also has ZPp = and as a reference spectrum uses 
a constant in f\ normalized in such a way as to have the same flux as 
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Table 1. Am for Johnson-Cousins and Stromgren filters using as red and 
blue spectra solar-mctallicity main-sequence Kurucz models with T c s of 
3 500 K and 50 000 K, respectively. 



Vega at the pivot wavelength of the Johnson V filter, 5493 A. We obtain 
a value of f\ = 3.5 8 2 • 10 ~ 9 erg s _1 cm~ 2 A -1 from the Vega spectrum of 
iBohlin fc Gilliland ( 200|). 



The ABMAG system also has ZPp = and as a reference spectrum uses 
a constant in f v normalized in such a way as to have the same flux as 
Vega at the pivot wavelength of the Johnson V filter, 5493 A. We obtain 
a value of f„ = 3.605 ■ 10~ 20 erg s _1 cm -2 Hz -1 from the Vega spectrum 
of lBohlin fc Gilliland (|2004h . 



• The Johnson-Cousins system uses Veg reference spectrum and has 

values of ZP P for P = U, B, V . . . that NEED TO BE MEASURED. 

There are some possible error sources that can yield biases when using 
synthetic photometry. The first one is using the wrong equation (^or[2I). Since 
Eq. n gives more weight to longer wavelengths, it generates brighter magnitudes 
for red objects than for blue ones compared to Eq. |2J This can be quantified 
e.g. by computing Am = (mp t i — m' P1 ) — (mp^ — m'p 2 ) for a red (1) and a 
blue (2) SED (Table since that quantity gives the maximum systematic error 
introduced in the analysis of a CMD. If one wants to convert between energy- 
integrating and photon-counting magnitudes it is useful to define a sensitivity 
curve .P'(A) = P(A)/A. Then, it is easy to show that mp> = m' P if ZPp/ = ZPp. 

A second possible source of errors is an incorrect /x r P f(A). For the case 
of Veg a, the most precise spectrum is the one obtained bv IBohlin Sz Gillilandl 
(|2004h using a combination of HST/STIS CCD spectroscopy and Kurucz models. 
The abs olute flux cal ibration has an accuracy of 4% in the FUV and 2% in the 
optical iBohlin 200(1} but, g iven that the photometric repeatability of STIS is 
0.2-0.4% (iB ohlin et al.ll2001l ). the relative flux calibration for colors derived from 
STIS spectra is expected to be better than 2% in the optical. 

Two other possible sources of error are an incorrect knowledge of P(X) and 
of ZPp. In the next two sections we analyze the cases of the Tycho-2 B^Vr and 
Johnson UBV filter sets. 



2. Calibration of Tycho-2 BtVt photometry 

The Hipparcos mission (|ES Al H997i ) observed the full sky and yielded the 
Tycho catalog, which is reasonably complete do wn to V = 11.5. The original 
Tycho catalog was consequently reprocessed by Hog et al.l (|2000bh to produce 
the Tycho-2 catalog, which contains 2.5 million stars and is currently the most 
complete and accurate all-sky photometric survey in the optical. The Tycho-2 
catalog contains photometry in two optical bands, Bt and Vt, whose sensitivities 
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Figure 2. Comparison between photometric and spectrophotometric Tycho- 
2 Bt — Vt colors as a function of photometric Bt — Vt for the two samples. 
The error bars represent the photometric uncertainties and the horizontal line 
marks the proposed ZPb t -v t - 



were analyzed by iBesselll ({2000) and use Veg the refere nce spectrum. Th e 
careful processing of Tycho-2 photometry was described by iHcte et all (l2000aft . 
including the differe nt tests used to check for possible systematic errors. 

I have recently (|Maiz Apellaniz! I200 5) tested the sensitivity curves of the 
Tycho-2 photometry and calculated the values of ZPb t , ZPy T , and ZPb t -v t us- 
ing the spectra obtained with STIS for the Next Generation Spectral Library (or 
NGS L, http://l if shitz .ucdavis . edu/~mgregg/gregg/ngsl/ngsl .html and 
iGregg et al. 2004) and the spectrophotometric standards of Bohlin et al. (2001). 
The results for Bt — Vt are shown in Figs. [21 and El No general trend is observed 
as a function of color in Fig. [21 and the data are symmetrically distributed around 
a central value, which I take to be ZPb t -v t - I measured ZPb t -v t by calculating 
the weighted mean using l/cf} y T as weights and found it to be 0.020±0.001 
magnitudes. The histograms for the Bt — Vt data, both in absolute and relative 
(corrected for ZPb t -v t and dividing each point by its photometric uncertainty) 
terms, appear in Fig. [31 The second histogram has a median of 1.2 • 10~ 5 and a 
standard deviation of 1.04 and the distribution is very well approximated by a 
normalized Gaussian. All of the above implies that an accurate cross-calibration 
of colors vs. relative fluxes between Tycho-2 photometry and HST spectropho- 
tometry is possible in principle without having to invoke e.g. modifications in 
the Tycho filter sensitivities or the STIS calibration. Furthermore, given that 
the normalized histogram has a standard deviation only slightly larger than 1.0, 
the largest source of deviations from the expected value originates in the pho- 
tometry, not in the spectrophotometry. Since the mean photometric ob t -v t = 
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Figure 3. Histograms for the comparison between photometric and spec- 
trophotometric Bt — Vt colors for the NGSL sample, (left) Regular histogram, 
(right) Histogram for the data shifted by the proposed ZPs T _v r and normal- 
ized by the individual uncertainties. A Gaussian distribution with [i = and 
a = 1 is overplotted for comparison. 



0.025 magnitudes, the accuracy of the spectrophotometrically-derived Tycho-2 
colors must be better than 1%, which agrees with the published value for the 
STIS photometric repeatability iBohlin et aD (|200ll ). A similar analysis for the 
individual magnitudes Bt and Vt yields ZP^ T = 0.078±0.009 magnitudes and 
ZPy T = 0.058±0.009. 

These results imply that the Tycho-2 photometry has well-characterized 
sensitivity curves and that the zero points are accurate and stable, both in 
absolute and relative terms. Therefore, it should be possible to produce accurate 
comparisons between Tycho-2 photometry and synthetic photometry generated 
from SED models. 



3. Calibration of Johnson UBV photometry 

The I Johnson! (|l966l ) UBV system is the most commonly used photometric 
system in the optical range. It was originally denned using photomultipliers but 
it was later adapted for CCDs. Despite its extensive use, its calibration for a 
comparison with synthetic photometry has a number of problems: 

• As opposed to the Tycho-2 case, the Johnson photometry available in the 
literature has been obtained by different observers using different detec- 
tors, telescopes, and observing sites, and applying different reduction pro- 
cedures. Under those circumstances, an artificial scatter of the measured 
magnitudes and colors is inevitable. 

• Ground-based photometry should not be as stable as space-based pho- 
tometry due to the larger instability of the observing conditions and the 
difficulties in correcting for atmospheric extinction, a problem that is ex- 
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acerbated for broad-band photometry with respect to the intermediate- or 
narrow-band case due to differential effects with wavelength. This point 
is especially important for the U band, whose short-wavelength limit is 
determined mostly by the atmosphere. 

• It is not clear that all published synthetic photometry works have i ncluded 
the correct distinction between Eqs. ^ and |2] (see e.g. ISTScllH998l ). 

A consequence of the above problems can be s een in the sens itivity curves 
for UBV published bv iBuser fc Kurucd (Il978ft and iBesselll (Il990ft . First, both 



of those articles are unable to find a unique B sensitivity curve that is capable 
of generating the observed U — B and B — V colors. Instead, they resort to two 
different definitions of B, one for each color, a result that is clearly unphysical. 
Second, the U sensitivity curves of those two articles are quite different: for 
A-type stars, the effect on the derived synthetic magnitudes is rather small, 
but for early- type O-stars the difference in the synthetic U — B is almost 0.1 
magnitudes. 

In order to reduce the effect of the problems above, I am working on a 
recalibration of Johnson UBV photom etry following an approach similar to the 
one followed in iMaiz Apellaniz! (|2005ft for Tycho-2 data. I have collected the 
Johnson p hotometry for the non-v ariable stars in the NGSL using the Lausanne 
database ( Mermilliod et all Il997t ) and selected those stars with at least four 



different observations in U — B or B— V in order to account for different observing 
conditions. For the absolute calibr ation of the photometry, I have used the value 
of ZPy = 0.026±0.008 obtained bv lBohlin fc Gillilandl (|2n04T ). Here I report the 
preliminary results. 

• The measured B — V colors agree well with the synthetic ones derived from 
the sp ect rophotometry and the sensitivity curves of either IBuser h Kurucd 
( 1978T I or IBesselll (Il990h. ZPr_v is very close to zero, in agreement with 



the result of iBessell et al.l ((1998). 



No combination of the IBuser fc Kur^Icl (|l978f ) and IBesselll (|1990T ) U and 



B sensitivities is capable of reproducing the observed photometry. 

I have used a ^-minimization algorithm to generate a new U sensitivity 
curve which, when combined with the B curve used to calculate B — V 
colors, is capable of reproducing the observed U — B colors. Therefore, it 
is possible to eliminate the discrepancies between the measured data and 
the synthetic photometry and, at the same time, get rid of the need of a 
double definition for B. 



4. CHORIZOS: maximizing information content 

As previously mentioned, astronomers usually attempt to derive the N pa- 
rameters of individual stars (temperature, luminosity, metallicity, extinction. . . ) 
or of stellar populations (age, metallicity, extinction. . . ) by comparing the ob- 
served SCMDs with the synthetic magnitudes derived from SEDs. That is, they 
transform their model atmosphere outputs from the theoretical (or parameter) 
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Figure 4. U — B vs. B — V color-color plot for Z = 0.0, main-sequence 
Kurucz atmospheric models. The line with circles indicates the location of 
the unrcdenned values as a function of temperature, starting at T = 50 000 
K, with the circles marking those points where the temperature is a multiple 
of 5 000 K . The rest of th e lines indicate the colors as a function of reddening 
using the Cardel li et alJ l)1989f) law with -R5495 = 3.1, with the color code 
being used to differentiate among temperature ranges which are relevant to 
determine the number of possible temperature + reddening solutions for a 
given (U — B) + [B — V) combination. Nine measured stars are marked, 
each one of them with ajj = ob = &V = 0.026. Note that the reddening 
lines are not straight and that they are not parallel to each other because no 
Q-parameter approximation is used. 
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Figure 5. Likelihood contour plots produced by CHORIZOS for the nine 
stars shown in Fig. 0] The x axis corresponds to the temperature in K and 
the y axis to the reddening £(4405 - 5495). 
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N — dimensional space to an observational M-dimensional space with M = 2 
(one color and one magnitude) and compare the results with the data there. 
Such a strategy has a fundamental problems: if N > 2, degeneracies are likely 
to exist, so a single solution cannot be obtained (this can happen even if N = 2). 
Also, it is not straightforward to include information from an arbitrary number 
of additional colors: this is usually done by working on individual planes (e.g. 
projections onto two-dimensional observational space or color-color diagrams) 
but in doing so one loses the full multidimensional information. It is also com- 
mon in such cases to use linearizing approximations, such as Q-ratio extinction 
corrections or filter transformations, hence introducing additional systematic 

errors (see Fig. QJ. 

An alternative is to use a code such as CHORIZOS ( Mai z"^pellanizll2004h . 



CHORIZOS first calculates the theoretical colors for the full ^-dimensional 
parameter space and then evaluates the likelihood in each point of the grid for 
the observed photometry. The moments of the resulting likelihood are then 
calculated to derive the range of possible parameters. Furthermore, since the 
code does not operate by simply finding a likelihood maximum (or \ 2 minimum) 
it allows for the detection of multiple solutions. A simple example for the well- 
known case of finding the temperature and extinction (JV = 2) for stars of 
known metallicity and gravity and with a standard extinction law from Johnson 
U — B, B — V colors (M = 2) is shown in Figs. Q]and|5J Such an example has 
degeneracies, as it can be seen in the two possible solutions for e.g. star 2 and 
in the three possible solutions for e.g. star 3. Such degeneracies can be lifted by 
using additional colors (Fig. EJ), which CHORIZOS can easily handle. 

CHORIZOS is available to the astronomical community and can be down- 
loaded from http : //www. stsci . edu/~ jmaiz A number of upgr ades to the code 
have been added since the publication of the original article ( Maiz Apellanid 



2004) and more are planned for the near future. We detail the current and 
future capabilities here: 

CURRENT FUTURE 

SEDs: Kurucz, Lejeune, TLUSTY Arbitrary user-defined SED models 

(stars); Starburst99 (clusters). for stars, clusters, galaxies. . . 

Parameters: Up to four. Two SED- Five, user-defined and in any corn- 
intrinsic (T e ff + log<7, age) and two ex- bination (metallicity, redshift, IMF 
trinsic (reddening + extinction law). slope, upper mass limit. . . ). 

Parameter control: Full range, re- Use of Bayesian priors 
stricted range, fixed. Grid size ad- 
justable and extendable. 

Filters: 80 preinstalled (Johnson- Future HST instruments?, user- 
Cousins, Strdmgren, Tycho-2, SDSS, defined. 
2MASS, HST). 

Spectrophotometry: No. Yes. 
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Figure 6. (left) V - I vs. B - V vs. U - B 3-color plot for Z = 0.0, 
main-sequence Kurucz atmospheric models. The range plotted covers T = 
3 500 - 50 000 and £( 4 405- 5495) = 0.0 - 2.0 and the extinction law used is 
that of iCardelli et al.l l)1989|) with -R5495 =3.1. The color surface marks the 
location in 3-color space while the black one is the projection onto the (U — B)- 
(B — V) plane (see Fig. 0}. (right) Basic topology for an M = 3, N = 2 case 
such as the one on the left panel. Given that M > N, the measured colors 
(blue circle) always lies outside the solution surface (red grid). However, 
including the uncertainty ellipsoid (blue grid), yields an intersection surface 
(solid yellow) of likely solutions. 



Eventually, CHORIZOS will become a full Bayesian code that will be able 
to handle any arbitrary SED family or filter set. 

As an example of the possible applications of CHORIZOS to the study 
of stellar p o pulat ions we include here a summary of the study performed by 
lArias et al.l (|2005l ). Those authors analyzed UBVIJHK S photometry for six 
stars in the Galactic H'n region M8. CHORIZOS was executed using solar- 
metallicity, main-sequence Kurucz stellar models with three free parameters: 
temperature, reddening, and extinction law. CHORIZOS parameterizes red- 
dening by i?(4405 — 5495), the monochromatic equivalent of E(B — V) and 
extinct ion law by R^a^, the monochromatic equivalent to Ry. The extinction 
laws of lCardelli et al.l (1989) were used. 

The likelihood plots for SCB 325, one of the six stars is shown in Fig. [7J 
CHORIZOS produces as graphical output the projection of the likelihood into 
each of the possible combinations of two free parameters (in this case, there are 
three free parameters, yielding three possible combinations). As it can be seen, 
the output is a well-defined peak with only a slight asymmetry, so the values 
of T c g, i?(4405 — 5495), and -R5495 and their uncertainties can be derived with 
confidence. 

Another of the CHORIZOS outputs is shown in Fig. [SI for the case of Her- 
schel 36, another star in the M8 sample. The left plot shows the best-fit spectrum 
when all seven photometric bands (six colors) are used. Such a best fit has a re- 
duced x 2 °f 21, indicating that the observed photometry is not compatible with 
the SED family and/or the parameter range used. A look at the plot immedi- 
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Figure 7. Likelihood contour plots produced by CHORIZOS for one of the 
stars in the M8 sample. Each plot corresponds to the projection onto one of 
the 3 orthogonal planes defined using the three coordinates T e g, £(4405 — 
5495), and #5495: T cff - £(4405 - 5495) (upper left), T cS - i? 549 5 (lower left), 
and £(4405 - 5495) - i? 5 495 (lower right). 
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Figure 8. Measured photometry and best-fit spectrum for Herschel 36 pro- 
duced by CHORIZOS. The left plot shows the fit using all seven photometric 
bands (UBVIJHK S ) while the right plot uses only the first five (UBVIJ). 
The corresponding values of the reduced minimum y 2 are 21 and 1.7, re- 
spectively. For the photometry, the vertical error bars are used to show the 
photometric uncertainties while the horizontal ones show the approximate 
extent of each filter in wavelength. 

ately suggests the origin of the problem: the UBVI photometry is apparently 
well fitted but the NIR data is not. As it turns out, Herschel 36 has a large IR 
excess, so a Kurucz model is not a good approximation for its SED in the NIR: 
CHORIZOS detects such circumstance by yielding a reduced minimum x 2 ^> 1. 
If on the other hand, CHORIZOS is run excluding the HK S photometry (right 
plot of Fig. |SJ), the result is very different: the reduced x 2 f° r the best fit is 
close to 1, as it is readily apparent from the concordance between the spectrum 
and the photometry (note that the two rightmost photometric points are not 
included in the fit). 

Once the photometry of each star has been processed by CHORIZOS and 
the parameters (T cj j, £'(4405 — 5495), -R5495) and their uncertainties have been 
evaluated, it is possible to construct the theoretical HR diagram (Fig. IHJ). Note 
that the uncertainty ellipses are inclined with respect to the two coordinates. 
The reason is the strong correlation between T e g and Mbd induced by the de- 
pendence of the bolometric correction on temperature. The HR diagram shows 
that the photometry of the four earliest-type stars in the sample is compatible 
with a distance modulus of 10.5. The two late B stars, however, must be far- 
ther away, since their uncertainty ellipses fall clearly below the expected main 
sequence at that distance. 



5. Conclusions 

I have shown the importance of the accurate calibration of photometric sys- 
tems in order to produce meaningful comparisons between the observed colors 
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Figure 9. Theoretical ( T e s vs. Mboi + 51ogri — 5) HR diagram for the 
six stars in the M8 sample produced by CHORIZOS. Ellipses indicate the 
68% likelihood contour for each star. The line shows the main sequence for 
51ogri — 5 = 10.5. As seen in the plot, such a distance is compatible with the 
CHORIZOS output for the four earliest stars but not for the latest two. 
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and m agnitudes and model SEDs. The Tycho-2 calibratio n derived bv iBesseTi 
(2000) (for the sensitivity curves) and lMaiz Apellanizl ( 20051 ) (for the zero points) 



is precise enough for that task. On the other hand, the available calibrations 
for Johnson UBV photometry yield relatively large systematic errors, which has 
prompted me to develop a new, more precise calibration. I have also shown how 
the use of multicolor photometry has significant advantages over the standard 
single-color + magnitude diagrams, among them the elimination of degenera- 
cies, the inclusion of multiple parameters, the avoidance of linearizing approxi- 
mations, and the possibility of a precise treatment of errors. 
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